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ABSTRACT 

An  analytical  formulation  is  presented  for  normal  incident 
plane  wave  scattering  from  gratings  of  compliant  tubes  embedded 
in  a  viscoelastic  layer  immersed  in  fluid.  Approximations  to 
boundary  conditions  at  the  interfaces  between  the  interstice  and 
compliant  elements*  are  modified  by  treating  the  interstice  as  a 
viscoelastic  plate  which  has  the  material  properties  of  the 
polymer  matrix  material.  The  compliant  tubes  are  modeled  as  two 
structurally  and  acoustically  coupled  plates  with  end  conditions 
which  correspond  to  those  for  the  vibrational  modes  of  a  highly 
eccentric  long  cylindrical  shell.  Conditions  on  normal 
displacement,  normal  stress,  and  tangential  displacement  relate 
the  compliant  plates  and  viscoelastic  plate  boundaries  to  the 
surrounding  viscoelastic  layers.  Comparisons  of  the  insertion 
loss  for  gratings  in  fluid  and  gratings  in  viscoelastic  layers  in 
fluid  are  shown.  Limitations  of  the  analysis  are  discussed  and 
calculations  will  be  compared  with  experimental  data. 

This  paper  was  presented  to  the  106th  meeting  of  the 
Acoustical  Society  of  America,  held  in  San  Diego,  CA.  during 
November  of  1983. 


ADMINISTRATIVE  INFORMATION 

This  work  was  done  under  NUSC  project  A70204  "Multiple 
scattering  from  shells,"  Principal  Investigator  Dr.  R.  P. 
Radlinski  and  Navy  project  No.  ZN00001  ,  program  element  6 1 1 5  2  N , 
Capt.  Z.  L.  Newcomb,  Naval  Material  Command  (Code  05B).  The 
authors  of  this  memorandum  are  located  at  the  Naval  Underwater 
Systems  Center,  New  London,  CT  06320. 
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I .  Introducti on 

In  previous  presentations,  the  interaction  between  multiple 
gratings  of  compliant  tubes  in  fluid  has  been  discussed  for 
several  different  grating  configurations.  In  this  work,  the 
transmission  of  plane  waves  through  a  single  grating  of  the 
compliant  tubes  embedded  in  a  finite  thickness  viscoelastic 
material,  surrounded  by  fluid,  will  be  investigated  through 
extensions  of  the  fluid  model.  Plane  strain  viscoelastic  theory 
is  used  to  describe  the  elastomeric  layer.  The  structural  modes 
of  the  compliant  tubes  are  treated  as  two  coupled  plates  with 
vibrations  corresponding  to  those  of  an  highly  eccentric  oval 
shell.  Differences  between  modeling  the  interstice  between  the 
tubes  as  an  elastic  layer  or  a  thin  elastic  plate  are  shown  to  be 
material  dependent.  In  both  cases,  the  layer  thickness 
transmissions  resonances  resulting  from  shear  wave  generation  by 
the  compliant  tube  inclusions  are  found  to  be  significant  in 
elastic  layers.  The  introduction  of  damping  into  the  elastomer 
mitigates  the  transmission  peaks  due  to  the  thickness  resonances. 
Comparisons  of  theoretical  calculations  with  experimental 
measurements  of  insertion  loss  for  panels  with  low  and  medium 
stiffness  viscoelastic  material  are  found  to  be  in  reasonable 
agreement . 

II.  THEORY 

A  schematic  representation  of  the  compliant  tube  grating 
embedded  in  a  viscoelastic  layer  is  shown  in  figure  1.  A 
normally  incident  plane  wave  that  transverses  region  1  strikes 
the  elastic  layer  at  the  boundary  of  region  2  and  a  part  of  this 
wave  is  reflected. 

The  portion  of  the  incident  wave  which  enters  the  region  2 
undergoes  multiple  reflections,  at  both  the  tube  and  elastic- 
fluid  boundaries  and  is  either  absorbed  or  finally  emitted  from 
the  layer  to  either  region  1  or  region  5.  Because  of  a  net 
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volume  displacement  generated  near  the  symmetric  bending 
resonance  frequencies  of  the  plates,  the  tubes  perform  as 
acoustically  soft  inclusions.  Shear  waves  are  generated  even  at 
normal  incidence  through  the  bending  and  rigid  body  translational 
motion  of  the  tubes.  Because  of  the  symmetry  in  this  grating 
arrangement,  the  solution  reduces  to  satisfying  the  boundary 
conditions  between  any  two  planes  of  symmetry.  Since  the 
solution  in  the  intertices  between  tubes  differs  from  the  rest  of 
the  layer,  the  encapsultant  layer  is  divided  into  three  regions. 

The  boundary  condition  in  the  fluid  at  the  symmetry  planes 
is  that  the  transverse  velocity  is  zero.  In  the  elastomeric 
material,  the  additional  boundary  condition  is  zero  transverse 
shear  at  the  symmetry  planes. 

The  total  pressure  In  the  fluid  above  the  elastomer  (region 
1  in  figure  1)  is  the  sum  of  the  incident  plane  and  an  infinite 
series  of  reflected  waves  as  given  by  the  first  expression  in 
figure  2.  Rn  represents  an  unknown  scattering  coefficient.  The 
wavenumber  of  the  nu  reflected  wave,  namely  kn  ,  is  real  when  kQ, 
the  wavenumber  in  the  fluid,  is  greater  than  the  grating 
parameter  an  In  this  case  the  reflected  wave  propagates  away 
from  the  grating  into  the  fluid.  For  kQ  less  than  an>  kn  is 
imaginary  and  an  evanescent  wave  propagates  nearly  parallel  to 
the  grating  and  decays  exponentially  away  from  the  grating.  In 
the  assumed  linear  viscoelectic  layer,  the  scalar  and  vector 

displacement  potentials  for  the  dilatational  (0.)  and  shear  waves 

t  h  ^ 

(A.)  in  the  j  elastomer  region  are  given  as  infinite  series  of 

J 

inhomogeneous  waves  traveling  in  both  directions.  The 

dilatational  and  shear  wavenumbers  k  and  k  respectively  are 

P 

assumed  complex.  The  constant  Cj  is  determined  by  the  region 
geometry.  A  j, J' ) ,  ,  ci^  ,  and  Dp  ^ ,  are  scattering 

coefficients  determined  by  satisfying  the  boundary  conditions  at 
the  region  interfaces.  Later  the  interstice  (region  3)  will  be 
modeled  as  a  viscoelastic  thin  plate  and  compared  with  the 
interstice  modeled  as  a  viscoelastic  layer. 
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Implicit  in  the  solutions  for  the  viscoelastic  interstice 

layer  between  tubes  is  the  assumption  that  at  both  the  symmetry 

plane  bisecting  the  interstice  between  the  tubes  and  at  the 

tube-elastomer  abuttment,  u  and  o  are  zero.  The  assumption 

j  xy 

that  a  is  zero  at  the  tube  abuttment  implies  that  the  shearing 
xy 

stress  caused  by  the  motion  of  the  symmetric  clamped  end  modes  is 
negligible.  If  the  tubes  are  thin  enough  this  assumption  is 
reasonable.  Because  of  symmetry,  the  condition  on  u^  is 
explicitly  true.  The  guided  end  conditions  are  matched  by  the 
assumed  sol ut i on  . 

Of  particular  interest  in  this  presentation  is  the 
transmitted  pressure  in  region  5  of  figure  1.  For  grating 
spacings  less  than  a  wavelength  in  the  fluid,  the  transmitted 
pressure  in  the  field  is  determined  by  T  ,  which  Is  the 
scattering  coefficient  of  the  single  transmitted  propagating 
wave . 


The  plates  are  assumed  to  obey  tThe  thin  plate  theory.  The 
totally  symmetric  modes  of  the  compliant  tube  which  result  in  a 
net  volume  displacement  are  approximated  by  clamped  end 
conditions,  and  are  shown  in  part  I  of  figure  3.  The  asymmetric 
modes  shown  in  part  III  include  rigid  body  motion  and  are 
described  by  guided  end  conditions  or  zero  slope  and  zero 
transverse  shear.  Complete  descriptions  of  these  wave  functions 
are  described  in  references  2  and  3. 

When  the  interstice  is  modeled  as  a  viscoelastic  plate,  the 
boundary  conditions  are  chosen  to  match  those  for  the  tube  model 
shapes.  However,  as  seen  in  part  II  of  figure  3,  the  motion  of 
the  viscoelastic  plate  is  different  than  for  the  two  connected 
plates  in  the  sense  that  the  motion  is  asymmetric  with  respect  to 
the  major  axis  of  the  plate.  The  guided  end  conditions  result  in 
Identical  mode  shapes  of  the  two  coupled  plates  of  the  tube  and 
the  interstice  modeled  as  a  thin  viscoelastic  plate.  The 
limitation  of  modeling  the  interstice  as  a  thin  viscoelastic 
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plate  Is  the  small  range  of  viscoelastic  materials  for  which  the 

4 

thin  plate  theory  is  valid.  From  Junger  and  Feit  ,  the  frequency 
range  that  the  classical  thin  plate  theory  is  adequate  is 
restricted  to: 


<  0.1  .  (1) 


Here  c$  is  the  shear  wave  speed  of  the  material  and  h  is  the 
interstice  thickness.  For  frequencies  larger  than  the  above 
restriction,  rotational  inertia  and  transverse  shear  must  be 
included.  At  even  higher  frequencies,  compress i ona 1  wave  effects 
must  also  be  included  in  the  modeling  of  the  intertice.'  Damping 
effects  have  been  incorporated  into  the  plate  theory  by  assuming 
a  complex  Young's  modulus. 

The  mathematical  expressions  for  the  wave  functions  which 

describe  clamped  and  guided  end  conditions  are  shown  in  figure 

4.  Both  expressions  contain  even  terms  because  only  modes 

symmetric  with  respect  to  the  center  (x-axis)  of  the  tube  are 

excited  by  the  normally  incident  plane  wave.  The  resonance 

frequency  is  proportional  to  the  square  of  the  wavenumber.  The 

first  non-zero,  in-air  frequency  of  the  guided  modes  occurs 

approximately  1.8  times  the  value  for  the  first  totally  symmetric 

compliant  mode.  The  non-zero  frequency,  guided  end  modes  have 

been  shown  to  degrade  insertion  loss  performance  for  closely 

2 

space  multiple  gratings  . 

The  boundary  conditions  to  be  satisfied  at  the  various 
region  interfaces  are  displayed  in  figure  5.  The  fluid- 
elastomer  and  elastic-elastic  conditions  are  those  common  for 
plain  strain  problems.  The  continuity  of  tangential  displacement 
at  the  pi ate-el astometer  boundary  is  derived  by  assuming 
i nextenti onal  motion  of  the  plate  such  that 


tan 


h  dW 

7  3y 
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III.  NUMERICAL  EXAMPLES 

In  the  first  example,  the  encapsulant  is  treated  as  a 
elastic  material  with  zero  damping  for  two  values  of  the  shear 
modulus.  The  dilatational  modulus  and  density  of  the  elastic 
material  is  similar  to  that  of  water.  The  insertion  loss  defined 
as  20  log1Q(P^/Pt)  a  measure  tlle  energy  transmitted  through 
the  baffle.  A  large  insertion  loss  implies  that  the  transmitted 
energy  is  small.  The  frequency  axis  has  been  normalized  to  f ^ 
which  is  the  first  symmetric  flexural  resonance  of  a  single  tube 
in  air.  The  same  normalization  will  be  used  throughout  the 
remainder  of  the  paper.  The  thickness  of  the  interstice  layer  is 
about  0.013  Xj  where  is  the  acoustic  wavelength  at  f ^ .  The 
elastic  separation  between  the  tubes  is  about  the  one  tube  width. 
Using  the  equation  (1)  for  the  validity  of  the  thin  plate  theory 

Q  O 

that  indicates  that  for  the  higher  shear  modulus  (10  nt/m  )  the 
classical  plate  theory  should  be  a  good  approximation  to  3^. 

As  seen  in  the  right  hand  side  of  the  figure  6,  good 

agreement  is  found  between  plate  theory  and  elastic  layer 

calculations.  Although  not  shown,  the  agreement  is  consistent 

6  2 

down  to  0.25fj.  For  the  low  shear  modulus  material  (10  nt/m  ) 
classical  plate  theory  is  valid  only  below  0.25^  and  poor 
agreement  is  found  between  the  Interstice  (region  3)  modeled  as  a 
plate  versus  the  intertice  modeled  as  an  elastic  layer  with 
approximated  boundary  conditions  at  the  tube  elastomer  abuttment. 
The  conclusion  is  that  only  for  very  sti f f . el astomer s  will  the 
plate  model  be  valid  for  tubes  with  significant  separation 
distances  d.  Studies  with  closely  packed  arrays  indicated  the 
modeling  of  the  interstice  in  these  cases  is  not  sensitive  in 
the  sense  that  the  viscoelastic  separation  acts  simply  as  a  mass. 

In  the  second  example  showed  in  figure  7,  the  effects  of  the 
encapsulant  as  a  function  of  shear  modulus  are  shown  for 
densely  packed  compliant  tubes.  The  insertion  loss  for  an 
identical  unencapsul anted  grating  is  shown  by  the  solid  curve. 
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The  layer  model  is  used  for  the  interstice  but  as  mentioned 
above,  the  interstice  is  insensitive  to  the  two  modeling 
techniques  for  densely  packed  tubes.  The  maxima  in  the  insertion 
loss  at  about  1 . 2 f ^  and  6.0f^  are  due  to  the  first  two  totally 
symmetric  clamped  tube  modes.  The  exact  location  of  the  maxima 
are  determined  by  array  fluid  interactions.  When  the  tubes  are 
encapsulated  in  the  low  stiffness  elastomer,  the  location  of  the 
maxima  do  not  change  but  the  depth  of  the  null  decreases.  Part 
of  this  depth  of  the  null  difference  is  due  to  the  damping  factor 

O 

(5=0.2).  As  the  stiffness  is  increased  to  u  =  10  nt/m  and  the 
damping  factor  increases  to  0.5,  the  character  of  the  insertion 
loss  changes  dramatically.  The  additional  stiffness  on  the  tube 
increases  the  effective  resonance  frequency  of  the  structural 
modes.  The  stiffness  also  decreases  the  vibrational  motion  of 
the  layer  so  that  more  energy  tends  to  be  transmitted  through  the 
1 ayer . 


The  effects  of  damping  are  shown  in  figure  8.  In  this 
example,  an  elastic  layer  is  compared  with  the  stiffer 
viscoelastic  layer  of  figure  7.  The  large  oscillations  in  the 
insertion  loss  are  due  to  thickness  resonances  of  the  elastic 
layer  from  the  resulting  shear  waves  generated  by  the  motion  of 
the  tube.  The  damping  in  the  viscoelastic  material  converts 
these  shear  waves  to  heat  and  thus  smooth  the  response  of  the 
insertion  loss  curves. 

The  effects  of  the  thickness  of  the  elastomer  are  shown  in 
figure  9  for  two  thicknesses.  The  thickness  of  the  original 
layer  (tQ)  is  or  0.56X$  where  Xg  is  the  shear  wavelength  at 

f 1 .  When  the  thickness  of  the  encapulating  layer  is  reduced,  two 
effects  control  the  changes  in  the  insertion  loss.  First,  the 
decreased  thickness  of  the  elastomer  reduces  the  effective 
stiffness  surrounding  the  tube  walls.  Second,  thinning  the  layer 
will  shift  the  layer  thickness  resonances  to  higher  frequencies. 
In  figure  9,  it  appears  that  the  decreased  stiffness  is  the 
dominate  effect  which  not  only  lowers  the  insertion  loss  maximum 
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of  the  first  symmet ri c  mode  to  1 . 5f  ^  but  also  allows  a  higher 
amplitude  of  vibration  such  that  more  energy  is  reflected  at  the 
resonance  maxima. 

IV.  COMPARISON  OF  THEORY  AND  MEASUREMENTS 

In  the  last  two  figures,  the  calculations  for  the  insertion 
loss  for  two  encapsulants  with  tubes  densely  packed  are  compared 
with  experimental  data.  The  interstice  is  modeled  as  a 
viscoelastic  layer.  The  first  comparison  is  shown  in  figure  10 
for  the  encapsulant  with  the  low  shear  modulus  material. 
Differences  between  experimental  and  the  theoretical  calculations 
can  be  attributed  to  the  nominal  values  of  the  shear  values  used 
for  the  calculation.  The  material  measurements  made  with  the  8&K 
Young's  modulus  apparatus  as  discussed  in  reference  5  have  large 
uncertainties  for  low  shear  materials.  The  castable  B.F. 

Goodrich  pc  rubber  material  is  relatively  temperature 
i nsensiti ve. 

The  second  comparison  between  calculations  and  theory  is 
shown  in  figure  11  for  a  closely  packed  tube  grating  encapsulated 
in  the  high  shear  modulus  polyurethane.  Again  reasonable 
agreement  is  found  over  the  entire  frequency  spectrum. 

V.  CONCLUSIONS 

The  parameters  which  affect  the  transmissivity  of  energy 
through  an  encapsulated  compliant  tube  grating  include  the 
dynamic  compliance  of  the  compliant  tubes,  the  stiffness  and 
damping  factor  of  the  encapsulant  material,  and  the  thickness  of 
the  layer.  Low  stiffness  materials  with  some  damping  to  smooth 
the  thickness  resonances  do  not  significantly  alter  the 
unencapsulated  grating  performance.  Higher  stiffness  materials 
increase  the  frequency  of  highest  dynamic  insertion  loss  but  also 
decrease  the  dynamic  compliance.  The  layer  thickness  resonances 
are  most  pronounced  in  purely  elastic  materials.  For  a  given 
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thickness,  the  layer  resonances  were  more  prominent  in  stiffer 
viscoelastic  materials  in  which  the  longer  wavelengths  make 
absorption  of  shear  waves  more  difficult. 

In  the  frequency  range  where  the  classical  thin  plate  theory 
is  valid,  good  agreement  occurs  between  layer  and  plate  modeling 
of  the  interstice.  For  closely  space  gratings,  modeling  of  the 
interstice  is  not  sensitive  to  either  the  plate  or  viscoelastic 
layer  model.  However,  it  is  noted  that  modeling  the  interstice 
as  a  fluid  but  retaining  region  2  and  region  4  as  viscoelastic 
layers  gave  significantly  different  results  than  the  two 
interstice  models  discussed. in  this  presentation. 
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